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1. Introduction 

The study of the excited baryon spectrum is an important research program for lattice QCD. 
In this report we present preliminary findings for the nucleon excited spectrum using the Nf = 2 
anisotropic Wilson action. We see evidence for a J p = f state. 

We construct operators which transform according to the irreducible representations (irreps) 
of the lattice rotation group (the octahedral group for cubic lattices) [[]], |". There are six double- 
valued irreps: G\ g , H g , G 2g , G\ u , H u , and G211, where the "g" subscript denotes positive parity 
(gerade) irreps and the "u" subscript denotes negative parity (ungerade) irreps. We identify the 
continuum limit spins by subducing the continuum rotation group to the octahedral group. The 
pattern of continuum spin states in each octahedral irrep is summarized in Table [|. 

To extract the excited spectrum, we construct a large number of operators in each symmetry 
channel. It is necessary to include operators with quarks displaced with respect to one another to 
capture radial and orbital excitations. After selecting a manageable set of optimized operators, we 
use the variational method, diagonalizing the matrix of correlation functions, to extract the excited 
spectrum. To distinguish between scattering states and resonances, it is necessary to include multi- 
hadron operators in the matrix of correlation functions. However, we use only three quark operators 
in this work. 

Table 1: The number of occurrences of each double valued irrep A = {Gi,H,G2} of the octahedral group 
for different values of continuum J up to y . 
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2. Nucleon operators 

We construct baryon operators which transform as irreps of the octahedral group as described 
in |^| . The basic building blocks are covariantly displaced smeared quark fields (D >j , iff (x)) Aaa 
with flavor A, color a, and spin a. The smeared quark field y is displaced p links in the j direction 
(j = 0, ±1,±2,±3). Both quark smearing and gauge link smearing are necessary to reduce noise 
and coupling to high energy states [S]. We use Gaussian quark smearing and stout link smearing 
[Q]. From the displaced single quark operators, we construct elemental operators 

KfrJjkW = eabc{b^y(x))A a a(pf Y(x)) Bbj} (D { k p) \ff(x)) Ccr , (2.1) 

where e a b c is the antisymmetric Levi-Civita symbol. To construct nucleon operators we project to 
j _ l 7 _ 1 . 

We use several different patterns of displacements for our three quark operators, as summa- 
rized in Table ||. Taking linear combinations of these elemental operators, we project to the irreps 
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Table 2: Patterns of quark displacements employed in elemental operators. The displacement indices indi- 
cate the direction of the gauge-covariant displacement for each quark; i. j, k = {0, ±1 , ±2, ±3}. The quarks 
all have the same displacement length. 



Operator type Displacement indices 
single-site 

«>- 

singly-displaced 

• ® • i = i = —k k^O 

doubly-displaced-I ' J ' ' 



i=j=k=0 
i = j = 0,k^0 



doubly-displaced-L 

triply-displaced-T 

triply-displaced-0 



« = 0, \j\? \k\, jk^O 
i = -j, \j\ + \k\, jk^O 
\i\ + \ j\ ^ ijk^O 



Qj= ^J1±1L ==. (2.3) 



of the octahedral group. This results in hundreds of nucleon operators in each channel, which are 
then "pruned" to manageable sets as described in [gj. 

To optimize the variational method, we would like a set of low noise, linearly independent 
operators. We first eliminate noisy operators based on the signal to noise ratio in diagonal elements 
of the correlation matrix. To test for linear independence of a set of operators, we use the condition 
number of the normalized correlation matrix on time slice t = a t : 

Qj(a t ) 
yJCii(a t )Cjj{a t ) 

For completely orthogonal operators, the condition number is 1 while for completely degenerate 
operators, the condition number diverges. We therefore seek a set of operators that minimizes the 
condition number. Pruning first within each operator type and then across all remaining operators, 
we produce a set of 16 optimal operators in each channel. The operators used in this work are the 
same as those in [S]. 

We can write the correlator of two nucleon operators in terms of three quark propagators 

x Q^ f (x,f|x ,0), (2.4) 

where 2^|^( x ) f l x 0i0) is a single quark propagator of flavor A from source site xo at time t = to 
sink site x at time t. The nucleon correlation matrix is 

„ _ i -j f jz,(uud) ^(uud) _ „(uud) _ „{uud) 

~ C apY C m\ U (a\a)(p\p)m + U (a\p)(fi\a)(r\f) U (a\f)(p\p)(r\a) (a|0)0|?)(r|8) 
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~ u (p\P)(rm(a\?) ~ (p\a)(r\PK<x\r) (r\?)(W)(m (r\PKP\f)(a\a) r { ' 

where the c' a g^ are the coefficients to project the operators to irreps of the octahedral group. The 
"/" superscript identifies each individual nucleon operator used to construct the matrix. 

3. Computational method 

3.1 Variational Method 

We use the variational method [|], ||] to extract the excited spectrum from the matrix of corre- 
lation functions, numerically solving the generalized eigenvalue problem 

Cf{t)vf{t,t ) = a { n A \t, to )C^\t )v { ;\t,to), (3.1) 

where n labels the eigenstates and A the double valued irrep of the octahedral group. To reduce 
instabilities in the eigenvectors due to degeneracies and numerical uncertainties, we solve the eigen- 
value problem on a single time slice t*. This time is selected to be as small as possible to minimize 
the noise, but large enough that the eigenvectors have stabilized. The correlator on all other time 
slices is then rotated to this fixed basis of eigenvectors. 

The diagonal elements of the rotated correlation matrix Qj(t) are related to the energies by 

Cf \t) ~ e - E '^ (l + ^fV^)) +Y j a n {t)e- E " {, - ,0 \ (3.2) 

no- 
where 8E is the difference between E\ and the next closest energy. The summation is a correction 
which results from the fact that we diagonalize only at t* . On other time slices we expect a small 
overlap between the basis vectors and all other energy eigenstates. However, the coefficients a n (t) 
should be negligible near t* . We can calculate the fixed eigenvector effective energy via 



E e / f (t)= In 



(3.3) 



3.2 Lattice Action 



We used 24 3 x 64 anisotropic lattices with the temporal lattice spacing a t three times finer 
than the spatial lattice spacing. Gauge configurations were generated using the Wilson action. The 
scale was set with the Sommer parameter. We analyzed an ensemble of 430 configurations with 
M n = 400 MeV. 

3.3 Fit Details 

Because baryon operators create a baryon and annihilate an antibaryon and we use antiperiodic 
temporal boundary conditions, if we create a baryon state at t = 0, we also create a state which 
propagates backwards in time from the opposite temporal end of the lattice. Because fermions and 
antifermions have opposite intrinsic parity, the antibaryon propagating backward in time has parity 
opposite to that of the fermion propagating forward in time. 

Since the signals decay exponentially, the backward-propagating signal will only be above 
the noise level for time slices above some threshold value of time. If the forward signal has an 
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effective energy plateau in a region below this threshold then we can extract the energy of the 
state by fitting the correlation function with an exponential decay through the plateau without any 
interference from the backward propagating signal. This was, in fact, the case for all channels 
except for G\ u . In this channel, the backward propagating signal was dominated by the G\ g ground 
state, the lowest energy state in the spectrum. For our lattice, the backward-propagating signal 
decayed slowly enough and the temporal extent was small enough (due to the anisotropy) that the 
G\ u effective energies had significant backward contamination. To extract the energy levels, we fit 
the correlation functions as a forward exponential decay plus a backward exponential decay which 
was constrained by the G\ g ground state. 

We performed fully correlated % 2 minimization fits. For each state a fit range was selected 
such that the x 2 was minimized, the quality factor Q was maximized, the corresponding effective 
energy plot plateaued in the fit range, and the fit parameters were stable under small variations in 
the fit range. As stated above, each channel except for the G\ u was modeled as 

ct\t)=Ae- E ^^\ (3.4) 

The G\ u correlation functions were fit simultaneously with the G\ g ground state, constraining the 
Gig ground state energy to be equal to the energy of the backward-propagating G\ u state: 

Cf '») = Ae- £ . <G "' ) ( ? - f »)+Be- £ o Glg V-0 

= De- £ o Gls V*) (3 . 5) 

We estimated the uncertainty in the energy through a jackknife analysis. We repeated the fit for 
each single elimination jackknife sample. The average energy and its jackknife error are reported. 

3.4 Filtering 

The presence of the backward-propagating state in the G\ u channel made the Liischer method, 
the diagonalization of the correlation matrix on each time step, unworkable because it led to large 
numerical instabilities in the eigenvectors. We tested a method based on filtering out the backward 
signal prior to diagonalization. In a time interval where the backward signal is simply the ground 
state of the opposite parity channel we model the unrotated correlation matrix as 

C \f\t) = ^A n e- E ^ +Be- E ° C ( T -'°K (3.6) 

n 

We define the filtered correlator as 

C%M) = Cf{t)-cf{ h ) + {l-e-^) t Cf{j) (3.7) 

j=t+i 

( e -Ett-to)_ e -E^ 1 - to )y (3g) 

where t\ is a time where the backward signal is, in fact, described by a single exponential. The 
filtered correlators consist of the forward signal plus a constant term. The diagonalization of the 
filtered correlators using the Liischer method produced stable eigenvectors, and the energy of the 
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states could be extracted by fitting the principal correlation functions to a single exponential decay 
with a constant term. This method did not produce any significant improvement over the results 
from the fixed eigenvector method, but we point out that the filtering is necessary to extract the 
excited spectrum using the Luscher method when there is significant backwards contamination in 
the effective energy. 

4. Results 

Four states were obtained in each channel. A plot of the spectrum is shown in Fig. |I[ In the 
positive parity channels, we can identify the G\ g ground state with the nucleon. Due to the high 
degree of degeneracy in the excited positive parity channel, it is difficult to identify the excited 
states. 

Results for the low energy negative-parity excited states allow for some interesting interpreta- 
tions. In the G\ u channel, we see the energy of the lowest energy states are close to the sum of the 
pion and nucleon mass. This raises the possibility that one of these states is, in fact, a pion-nucleon 
state. Even if this is the case, the lowest negative-parity nucleon state in G\ u is below the first 
excited positive-parity state in G\ g . This does not match the physical spectrum where the lowest 
negative parity resonance, the /V*(1535) is above the first excited positive parity state, the /V(1440). 

In the H u channel, we expect the lowest state to correspond to a spin f state. We interpret 
the lowest H u state as the /V(1520). In the G% u channel we see that the lowest-energy is degenerate 
with a partner H u state, with no state in the G\ u channel at the same energy. This is the signature of 
a spin | state. Two of the six components needed for a spin | state are occur in the G2 U channel 
and four in the H u channel. The Gi u and H u states corresponding to the | must be degenerate in 
the continuum limit. Moreover, there should not be a G\ u state that is degenerate with these two 
states as that would indicate the possible presence of a spin Z state. These conditions are satisfied 
near Ea t = 0.35 giving clear evidence for a | state. We interpret this to be the lowest | state in 
the physical spectrum, the /V(1675). 

A | state has not been seen in any earlier work. Our quenched QCD analysis had three 
degenerate states (within errors): G2 U , H tl and G\ u . That pattern had two possible interpretations: 
a single spin | state or an accidental degeneracy of a spin | state and a spin ^ state. 
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Figure 1: The nucleon spectrum obtained for each symmetry channel for 24 x 64 Nf — 2 lattice QCD data 
at m K = 400 MeV. Errors are indicated by the vertical size of the box. 
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